The Warped Disk of Centaurus A from a radius of 2 to 6500 pc 
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Abstract: 

We compile position and inclination angles for tilted ring fits to the warped dusty and gaseous disk 
spanning radius 1.8 to 6500 pc from recent observations. For radii exterior to a kpc, tilted ring 
orientations lie on an arc on a polar inclination versus position angle plot, suggesting that precession 
following a merger can account for the ring morphology. Three kinks in the ring orientations are seen 
on the polar plot, the one at radius of about 1.3 kpc we suspect corresponds to the location where 
self-gravity in the disk affects the ring precession rate. Another at a radius of about 600 pc may be 
associated with a gap in the gas distribution. A third kink is seen at a radius of 100 pc. A constant 
inclination tilted disk precessing about the jet axis may describe the disk between 100 and 20 pc but 
not interior to this. A model with disk orientation matching the molecular circumnuclear disk at 100 
pc that decays at smaller radii to an inner flat disk perpendicular to the jet may account for disk 
orientations within 100 pc. Neither model would account for the cusps or changes in disk orientation 
at 100 or 600 pc. 

Keywords: Write keywords here 



1 Introduction 

In its central regions, NGC 5128 exhibits a well recog- 
nized, optically-dark band of absorption across its nu- 
cleus. This dusty dis k was f i rst mo deled as a transient 
warped structure bv iTubb^ ||1980D . The kinematics of 
the ionized and molecular gas are well modeled by a 
warped disk composed of a series of inclin ed connected 
rings un dergoing nearly circular m otion ( Bland 1 9861: 
Bland e t al. 198?!: iNicholson et al.l E992: Guillen e taP 
19921 1. Recent work has better fit the morphology of 
the disk in th e infrared using Spit zer Space Telescope 
observations (|Quillen et al ] l2006bl ) (and improving on 
the mod el fit to disk seen i n absorption in the near- 
infrared; iQuillen et al.|[T993 ). with high resolution in- 
terferometric molecular C O observations and cover- 
ing the central arcminute (|Espada et al.l [20091 ). trac- 
ing the warm molecular hydrogen to within the sphere 
of influence of the black ho le using adaptive optics 
(|Neumaver et all |2007|. l2009ri. a nd with new HI ob- 
servations l Oosterloo et al.ll2009l 'l. 

We take the opportunity to compile tilted ring fits 
to the disk morphology. The result is a warp ge- 
ometry, described by tilted ring major axis position 
and inclination angles as a function of radius, that 
spans from radius 0.1" to 3 50" or using a di stance 
of 3.8 Mpc to Centaurus A (|Harris et al.l 120091 ') from 
1.8 pc to 6.5 kpc. Here the smallest scales are set 
by the 0.12" K band resolution of the adaptive optics 



observations (jNeumaver et al.] l2007l ) and the largest 
scal es by the e xtent of the tilted ring fit to the HI 
( Oosterloo et al. .2009 '). The factor of 3600 between the 
smallest and largest scales is large compared to other 
studied gaseous disks. This large scale of observa- 
tions allows us to search for orientation changes caused 
by different physical processes. In the outer regions 
the primary processes causing orientation changes are 
torques from the galactic potential and from pertur- 
bations caused by galaxy mergers. Interior to these 
torques may be caused by a va riety of proce sses in- 
cluding self-gravit y in t he disk ('Spark^ |l996|'l. radia - 

tion pressure (e.g.. lPri"ng lc 1996; Malo nev et al.ll996l), 

winds or pressure gradients in the ambient ISM I Quill en fc Bowed 
il99S : lOuillen 2001.). and the Bardeen-Petterson effect 
(e.g., ICaproni et aLll2007l ). 

2 Compiling Tilted Ring Fits 

Position and inclination angle for the tilted ring fit to 
the warm molecular hydrogen emissio n 1-0 S(l) H2 line 
at 2.1 22/jm were taken from Table 1 bv lNeumaver et all 
(|2007l ) and have a spatial resolution of 0.12", covering 
the central 3" x 3". The submillimeter SMA CO obser- 
vations have a non-circular beam with wi dest direction 
6" F WHM and cover a 1' field of view (jEspada et all 
120091'). Angles us ed are those used to make Figure 8a 
bv lEspada et al.l (|2009l ') and are those for the tilted 
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Figure 1: Position and inclination angle as a tunc 
tio n of radius comp iled from table s and figures 
by iNeumaver et all (|2007l l2009l): lEspada et al ' 



(|2009[) : iQuillen et alj "(|2006b( ): lOosterloo et a! 
20091 ). Data used to make these plots are listed 



in Table 1. 



ring fit lacking a bar perturbation. The angles from 
the fit to the 8/im morphology seen in Spitzer Space 
Telescope imaging with t he IRAC camera are those 
used to make Figure 7 bv lQuillen etll] l|2006bl ). The 
IRAC camera at 8/im has a resolution of 2.2". The 
angles from the tilt ed ring fit to the new HI obser- 
vations are by Oost erloo etal] HqoI). There are no 
angles listed within 60" from the HI fit because HI is 
depleted in the center of the galaxy. The fifth set of po- 
sition and inclination angles are from observations and 
an associated tilted ring fit that will be described in 
more detail bv INeumaver et al.1 (l2009l 'l. The measure- 
ments are from integral field spectroscopy with SIN- 
FONI and are taken in natural seeing under excellent 
conditions of 0.5", cov ering a 8" x 8" field of view (see 
ICappellari et al.l |200£ ) and using the same 1-0 S(l) 
H2 line discussed by Neumaver et al.l (120071 ). Tilted 
rings were fit to the velocity fie l d usin g the procedures 
described by INeumaver et al.l (|2007l ). The compiled 
position and inclination angles from these five studies 
as a function of log radius in pc are shown in Figure 



[T] These angles are also listed in T able 1. We note 
that t wo sets of these tilted ring fits ([Neumaver et al] 
l2009l : lOosterloo et alJ l2009t ) are preliminary and may 
be updated in future. 
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Figure 2: Polar plot of ring position and inclina- 
tion angles. Here points are plotted with radius on 
the plot given by the ring inclination and at an an- 
gle on the plot set by the ring major axis position 
angle. The origin (with zero inclination) is on the 
upper right. The same data are shown here as in 
Figure [TJ We note that the disk near the nucleus 
lies at quite a different orientation than the outer 
disk. Three kinks or cusps are identified on this 
plot. The first on the lower left and in both the 
HI and Spitzer points (turquoise circles and pink 
squares). The second on the upper left and pri- 
marily in the CO points (blue stars). The third 
between the CO and H2 fits on the lower right 
(blue stars and green X-s). 



The inclination angle is defined such that an edge- 
on disk is oriented with i = 90° . Inclinations above 90° 
refer to those with northern side closer to the viewer 
(and so masking the northern side of the galaxy) and 
those below 90° refer to the opposite. The tilted ring 
fit to the disk seen in extinction (e.g., iQuillen et al.l 
I1993I ') specify which side of the disk is closer to the 
viewer. The position angle of a ring is that on the sky 
(counter clockwise from north) of the blue shifted side 
of its major axis. 

There is consistency between the different tilted 
ring fits in the outer galaxy with the HI based kine- 
matic fit at lower inclination than that done from the 
Spitzer images. The fit to the Spitzer images was done 
without kinematic constraints. Because the model was 
created with a 3D data cube, it was computed faster 
when the cube was smaller. This likely biased the best 
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fit model to configurations more twisted and compact 
than the real outer disk. The HI kinematic fit should 
be considered more accurate than the model used to 
match the Spitzer images. 

The CO data appear to connect the outer and inner 
regions though we must keep in mind that the molec- 
ular gas distribution is dominated in the center by a 
100-200 pc radius circumnuclear disk and the gas dis- 
tribution may contain one or more gaps. The inner- 
most points of the CO ring fit arise from the circum- 
nuclear disk itself and so have well measured position 
angles. Consequently it is unlikely that the innermost 
position angles have been overestimated. 

Th e fit s to the H2 emi s sion b y iNeumaver et al.l 
(|2007l 'l and iNeumaver et al.l (|2009l ') overlap and have 
consistent position angles, but those from the smaller 
field of view data have about 20° higher inclinations. 
W e expect that the fits to high resolution observations 
bv INeumaver et al.l l|2007t) would be more accurate in 
the center, and those to observa t ions c overing a larger 
field of view bv Neu maver et ahl (|2009l ) would be more 
accurate at larger radii. However further analysis (and 
possibly deeper observations) are required to resolve 
the discrepancy. 

2.1 Outer disk 

In Figure [2] we show the same angles listed in Table 1 
and shown in Figure [1] but on a polar plot. Here the 
radius of points on the plot is set by the inclination 
angle and the angle of points on the plot set by the 
position angle of the tilted rings on the sky. We first 
discuss the disk exterior to 1 kpc. We see that the 
ring angles from the outer galaxy lie on an arc. Were 
we to plot the ring orientations on a sphere the arc on 
this plot could correspond to part of a circle centered 
about a particular orientation. The rings in the inner 
galaxy lie at quite different orientations compared to 
the rings in the outer galaxy. 

In Figure [3] we show the same polar plot as in Fig- 
ure[2]but with two separate ovals. Each oval shows vec- 
tors that are the same angular distance from a particu- 
lar orientation. The red oval (on the left) corresponds 
to all rings that have angular momentum axis 20 de- 
grees from a direction that is 105° inclined from the 
line of sight and at a position angle of 18° on the sky. 
This is similar to the axis used t o set the underlying 
galaxy symmetry axis assumed bv lOuillen et al.l (|l99i 
l2006m when generating precession models for the disk, 
though in neither model was the disk inclination held 
fixed with respect t o the galaxy axis. The model by 
iQuillen et al.l 1119931 ) has inclination angle with respect 
to the assumed underlying galaxy symmetry axis that 
increases with increasing radius whereas that used to 
fit the Spitzer data decreases with increasing radius. 

The arc in the polar plot from the outer disk can 
be accounted for with a merger model. Following a 
merger, gas accreted from another galaxy would have 
angular momentum that is not necessarily aligned with 
the symmetry axis of the underlying galaxy. In a coor- 
dinate system with respect to a galaxy symmetry axis 
(assuming the galaxy is approximately axisymmetric) 
each gas ring would maintain its inclination but would 
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Figure 3: Same as Figure [5] but with the addition 
of two ovals and a black spiral. The ovals show 
orientations lying at fixed angular distance from a 
particular orientation on the sky. The spiral shows 
a disk that approaches an orientation perpendicu- 
lar to the jet axis at small radius. The outer disk 
orientations lie on an arc which can be explained 
with precession following a merger. Precession oc- 
curs about a symmetry axis that is estimated from 
the center of the red oval (on the left). The gray 
oval (on the right) shows ring orientations that 
have angular momentum axis 18° from orientation 
with inclination 45° and position ang le the same 
as the radio jet (55°: lBurns et al.1ll98 3*). The black 
spiral shows a model with inclination and position 
angle the same as the circumnuclear disk at 100 
pc but that approaches orientation perpendicular 
to the jet axis at small radii. 



precess about the galaxy symmetry axis. The angular 
difference between the gas angular momentum and the 
galaxy symmetry axis sets the size of the oval whereas 
the galaxy symmetry axis orientation determines the 
center of the red oval (on the left) in Figure |3l 

A simple precession model will not give a perfect 
fit or prediction because mergers are not expected to 
leave a gas remnant restricted to a single plane (cor- 
responding to a very narrow distribution of angular 
momentum orientation), the underlying galaxy itself 
may be triaxial rather than axisymmetric and would 
not be static during the merger. The outer HI ring 
orientations (r > 1.6 kpc or 85") lie shifted but not 
too distant from ring orientations in the disk between 
0.8 -1.5 kpc (or 40 - 80"), suggesting that the angular 
momentum vector of gas stripped in the outer regions 
was not too different from that forming the middle re- 
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gions of the disk. 

Within about 1.3 kpc the position angle increases 
with increasing radius rather than decreases with ra- 
dius as it does exterior to 1.3 kpc (Figure[T^). This is 
also seen as a kink or cusp in the points on the lower 
left in the polar plot (Figure [2]| . A tilted ring pre- 
cesses at a rate proportional to the ellipticity of the 
gravitational potential and the angular rotation of a 
particle in a circular orbit. Rotation rates are faster in 
the center of the galaxy so a warped disk should be in- 
creasingly corrugated in the center. A reduction in the 
rate or reversal of direction of precession can cause a 
cusp in the orientation angles. This kink would either 
be caused by a variation in galaxy eccentrici ty, but 
could also be due to self-gravity in the disk l|Sparkel 
119961 '). The direction of precession in the outer galaxy 
and isophotal major axis at a PA of about 20° suggests 
that either the galaxy is prolate and the disk is oriented 
near in its symmetry plane or the galaxy is oblate and 
the disk is in a nearly polar orbit (lOuillen et al.|[l993l : 
ISparkeiri996l ). A sharp change in galaxy eccentricity is 
not seen in the 3fim isophotes suggesting that future 
warp models should reconsider the importance of self- 
gravity as well as updated mass models for the stellar 
mass distribution in accounting for the disk morphol- 
ogy- 

The outer disk at 8 and 24 fim and in HI emission 
is not symmetrical abo ut the origin (e.g., see Figure 2 
bv lQuillen et al.ll2006bl ) implying that the outer disk is 
somewhat lopsided. The degree of lopsidedness could 
be precisely measured from the HI kinematic observa- 
tions. By studying the evolution of this asymmetry 
it might be possible to derive constraints on the time 
since the merger and on the nature of the merging 
event itself. 

2.2 Inner Disk 

In the polar plot (Figure [2]| there are at least 3 loca- 
tions where there are abrupt changes in orientation. 
One change occurs at a radius of about 35" or 600 pc, 
the radius at which lOuillen et al.l (|2006bl ) failed to find 
a good model for the Spitzer images and so suggested 
that there might be a hole in the gas and dust distribu- 
tion. lEspada et al.l (|2009l ) found evidence for a gap in 
the gas distribution between a radius of 200 and 800 pc 
and also showed that the fits could be improved with a 
model that included non-circular motions the gas. In 
Figure [5] the radius of the gap in the gas distribution 
corresponds to the topmost plotted point where there 
is a cusp in the CO generated points. 

There is a third cusp on the polar angle plot in 
Figure [2] on the lower right corresponding to orien- 
tation angle of the circumnuclear molecular disk at 
a radius of 100-200 pc (6-12"). Orientatio ns based 
on th e near-infrared spectroscopy by Ncuma ver et al.l 
(|2009l ') approach that of the circumnuclear molecular 
disk, suggesting that there is cusp in the angular ori- 
entations at a radius of about 100 pc or 6". 

The CO circumnuclear molecular disk could be 
smoothly connected to the disk seen in warm molec- 
ular hydrogen within 100 pc. Both the circumnuclear 
molecular disk at radius 100-200 pc and the innermost 



rings at a radius of 2 pc are approximately perpen- 
dicular to the jet axis, though the innermost points 
have lower inclinations than the circumnuclear disk. 
The radio and X-ray jets have a pos ition angle of 55° 
(|Burns et al.lll983l : iKraft et al.ll2000D with the north- 
ern jet suspected to be closer to the observer. VLBI 
studies suggest that the jet axis lies between 50 and 
80° from our line of sight (Tingav et a,l, 1998) , however 
apparent motions suggest an angle between 20 and 40° 
(jHardcastle et all [2003 ). A disk perpendicular to the 
jet would have a major axis position angle of 145°. If 
the disk at the smallest measured radii is perpendicu- 
lar to the jet then the jet is about 40° from the line of 
sight. 

On Figure [2] the points bv iNeumaver et all (|2009l ') 
(green X-s) are at hig her inclination than those by 
INeumaver et al.l l|2007l ) (red plus signs). Both fits at 
r ~ 80 pc begin near the orientation of the circumnu- 
clear disk on the left in Figure[2]and Figure [3l With de- 
creas ing radius, the green points (X-s) (Ne umaver et al.] 

|2pa9f) move upwards (to lower inclination) a nd then to 

the r ight, whereas the red points (plus signs) (|Neumaver et al.] 
l2007l l move to the right and then back to the left at 
the smallest radii. 

Perhaps the circumnuclear molecular disk orien- 
tation and that seen interior to it can be connected 
with an arc corresponding to precession about an axis 
near the radio jet. Orientations 18° from an angu- 
lar momentum axis inclined 45° from the line of sight 
and at a position angle of 55° are shown as a gray 
oval (on the right) in Figur e [3l The gray oval passes 
near most of the points bv .Neumaver et alj (2009) but 
is inconsistent w i th th e orientation in the inner disk 
(|Neumaver et al.]|2007t) . The jet inclination a ngle is 
uncertain (see discussion bv iHardcastle et ahl 12003). 
but even if we changed the estimated jet inclination 
angle and so moved the center of the oval, it would 
not pass through all the data points interior to the 
circumnuclear molecular disk in Figure [3l part i cularl y 
those at radii within 10 pc bv INeumaver et al] l|2007t) . 

We have also plotted on Figure|3]a black spiral that 
has an orientation near the circumnuclear disk and ap- 
proaches an orientation perpendicular to the jet axis at 
a radius of 1 pc. For this model, the angular distance 
between the angular momentum of the circumnuclear 
disk at 100 pc and the jet axis is 25°. The angle 
between the angular momentum vector and jet axis 
smoothly drops to zero with a slope di/da = —0.05 
where i is inclination and a precession angle with re- 
spect to the jet axis, as the precession angle advances 
by 450°. This spiral was projected assuming that jet 
axis is 40° from the line of sight. This model passes 
near data points at radii interior to 10 pc, suggesting 
that simultaneous precession and damping to a plane 
perpendicular to the jet axis may be occurring in the 
disk at radii interior to 100 pc. 

The gravitational field of the black hole dominates 
that from the stellar component within a radius known 
as the sphere of influence, Vs = GMbh /o-^ ■ For a black 
hole mass of Mbh = 5.5 ± 3.0 X IO^Mq and velocity 
dispersion of about 150 km/s (jCappellari et al.ll2009l ) 
this radius is about r s ^ 10 pc, consequently the inner- 
most data points bv INeumaver et "ahl l|2007t ) lie within 
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the massive blaclt hole's sphere of influence. For com- 
parison, the gravitational radius of the black hole at 
the nucleus is Vg = 2GMbh/c^ = 5 x 10~® pc. Rel- 
ativistic frame dragging near rg for a spinning black 
hole leads to precession known as Lense-Thirring pre- 
cession and warping of a disk that is initially mis- 
aligned with the black hole spin axis. This is called the 
Bardeen-Petterson effect. Due to viscous processes, 
the disk could settle into the plane perpendicular to 
the black hole's spin axis as it accretes. The transi- 
tion or Bardee n-Petterson radius i s estimated to be 
of order W^Vg (jCaproni et al.l [20071 ) corresponding to 
0.1 pc fo r Centaurus A. This is well within that re- 
solved bv Ncumaver et al.l (|2007t) . A similar situation 
exists for NGC 4258 with only the inner radius of 
the resolved warp ed disk near the estim ated Bardeen- 
Petterson radius (|Caproni et al.l [20071 '). The discrep- 
ancy between the estimated Bardeen-Petterson radius 
and the observed location of the warp might be re- 
solved if there is a smooth transition between the outer 
disk and the Bardeen-Petterson radi us, as exhibited by 
simu lations of diffusive disks (e.g., iLodato &i Pringla 
120061 '). The warped disk of NGC 4258 has been inter- 
preted in terms of the Bardeen-Petterson effect even 
though the warped disk itse lf lies outside the est imated 
Bardeen-Petterson radius ([Caproni et al.l [20071 ). The 
rough correspondence between the black spiral and 
data points in the inner disk in Figure |3] suggest that a 
similar model might be applied to Centaurus A's disk 
within a radius of 10 pc. Such a model would be a bet- 
ter fit for a jet axis near 40° from the l ight of sight, con- 
sisten t with the range estimated by iHardcastle et al] 
(|2003l '). bu t 10° lower th an estimated from VLBI ob- 
servations l|Tingav et al.l[T99 8). 

Outside the sphere of influence of the black hole, 
Vs, disk settling or warp diffusion timescal es exceed 
the a ccretion timescale ( Steimari-Cameron fc DurisenI 
Il988l l. However inside the black hole's sphere of influ- 
ence the rotation is nearly Keplerian. The warp dif- 
fusion timescale is significa ntly shorter for Keplerian 
disks than for galactic disks ijPapaloizou fc Pringlell98^ 
The warp diffusion timescale for a Keplerian disk td^ ^ 
a(^)"^ (jPapaloizou fc Pringld Il983l ') where h/r is 
the disk aspect ratio and a is the unitless parameter 
used to characterize the disk viscosity. The Bardeen- 
Peterson radius (at which a disk would settle to the 
midplane defined by the black hole spin axis) is es- 
timated by equating td to the inverse of the Lense- 
Thirring precession rate. The ratio of the Lense-Thirring 
precession rate to the angular rotation rate is at most 
10~* at 1 pc from the black hole. The expression for 
td illustrates that the Bardeen-Petterson radius could 
only be located near 1 pc if the disk is extremely thin 
and viscous, so much so that it is uncomfortable to ap- 
ply such a model to the inner region of Cen A's disk. 
One could consider the possibility that the disk does 
settle to a plane perpendicular to the jet axis and it 
does so near and within Vs where the warp diffusion 
timescale is shorter than exterior to this radius. If this 
is the case then there are other torques acting on the 
disk that depend on orientation angle, for example due 
to radiation pressure or pressure gradients that cause 
precession faster than Lense-Thirring precession near 



A model with a smoothly changing orientation an- 
gle that decays to an axis perpendicular to the jet 
would not explain the orientations at radii between 
100 and 600 pc, where there may be a gap in the gas 
distribution. Neither would a Bardeen-Petterson effect 
model account for the cusps or changes in orientation 
at a radius of 100 pc and at 600 pc. One possibility 
is that a triaxial bulge induces non-circular motions 
in this re gion complicating t he measurement of ori- 
entations (jEspada et al.l [20091 ). Another possibility is 
that bursts of star formation or nuclear activity have 
blown out material from this region (e.g.. Ouille n et al.l 
[2006 a). In this case, the kinks in the orientation direc- 
tion may be related to past episodes of activity. This 
last possibility suggests that future work could explore 
a relation between orientation and episodic nuclear ac- 
tivity even at radii interior to 100 pc. Orientation 
changes related to epochs of nuclear activity would be 
predicted if either radiative or wind driven warp insta- 
bilities operate or if there are variable gradients in the 
ambient medium associated with jet propagation. 

3 Discussion 

Surprisingly a comparison of the tilted ring fits shows 
at least three kinks in the orientation angles as a func- 
tion of radius on a polar plot. The first ( at a radius 
abou t 1.3 kpc or 70") was previously known (jOuillen et al.l 
Il993l ) and could correspond to a variation in precession 
rate caused by self-gravity of the disk. The outer disk 
forms an arc on a polar plot of ring orientation angles 
implying that p r ecessi on following a merger, as pro- 
posed bv lTubbsl ||1980D . is a good explanation for the 
disk morphology. 

Between the hole in the dust distribution at a ra- 
dius of about 30" (or 600 pc) and the circumnuclear 
molecular disk there is also an abrupt change in an- 
gle. There is another change in orientation at a ra- 
dius of about 100 pc (or 6"). Precession about the 
jet axis could match the orientations between 100 and 
20 pc, but would not match the orientations interior 
to this. A model with disk angular momentum axis 
approaching the jet axis with decreasing radius could 
be promising. Models of the Bardeen-Petterson effect 
that are used to interpret NGC 4258's disk might un- 
der extreme conditions be applied to Centaurus A even 
though the resolved region (outside 1.8 pc) lies out- 
side the estimated Bardeen-Petterson radius. How- 
ever, the two inner cusps or changes in the disk orien- 
tation would not be predicted by a Bardeen-Petterson 
type of model. Future work could explore models relat- 
ing disk orientation to past episodes of nuclear activity. 
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Table 1: Inclination and Position Angles from 
Tilted Ring Fits 
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"Angles are from T able 1 by In eumaver et al.1 
(200 f). 

''Preliminary fits by Neumaver et al.l (l2009t). 

^Angles are from F igure 8a by lEspada et al.l 
(I2009D . 

''Angles are from Figure 7 by iQuillen et al.l 
(l2006bD. 

'^Preliminary fits bv I Poster loo et al] ( 2009f l. 



